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Polymerization of actin proteins into dynamic structures 
is essential to eukaryotic cell life. This has motivated a large 
body of in vitro experiments measuring polymerization ki- 
netics of individual filaments. Here we model these kinetics, 
accounting for all relevant steps revealed by experiment: 
polymerization, depolymerization, random ATP hydrolysis 
and release of phosphate (Pi). We relate filament growth 
rates to the dynamics of ATP-actin and ADP-Pi-actin caps 
which develop at filament ends. At the critical concentra- 
tion of the barbed end, c crit , we find a short ATP cap and 
a long fluctuation-stabilized ADP-Pi cap. We show that 
growth rates and the critical concentration at the barbed 
end are intimately related to cap structure and dynamics. 
Fluctuations in filament lengths are described by the length 
diffusion coefficient, D. Recently Fujiwara et al. [Nature 
Cell Biol. (2002) 4, 666] and Kuhn and Pollard [Biohys. J. 
(2005) 88, 1387] observed large length fluctuations slightly 
above c cri t, provoking speculation that growth may proceed 
by oligomeric rather than monomeric on-off events. For 
the single monomer growth process we find that D exhibits 
a pronounced peak below c crit , due to filaments alternat- 
ing between capped and uncapped states, a mild version 
of the dynamic instability of microtubules. Fluctuations 
just above c crit are enhanced but much smaller than those 
reported experimentally. Future measurements of D as a 
function of concentration can help identify the origin of the 
observed fluctuations. 

The tendency of actin protein to spontaneously polymer- 
ize into rapidly growing filaments is fundamental to the life 
of eukaryotic cells. Cell motility [1,2], cell division [3], and 
endocytosis [4] are examples of processes exploiting the dy- 
namic character of actin structures composed of filaments. The 
regulation of filament growth processes leads to well-defined 
structures and coordinated function. For example, in combi- 
nation with branching, capping, and depolymerizing proteins, 
actin self-assembles into controlled dynamic cross-linked net- 
works forming the dynamic core of lamellipodia [2]. 

These complex cellular actin-based systems exhibit multi- 
ple superposed mechanisms. This has inspired a large body 
of in vitro work aiming to unravel these mechanisms and pin 
down rate constants for the constituent processes in purified 
systems [5], An important class of experiments entails mea- 
suring growth rate at one end by microscopy [6-9] or by bulk 
spectroscopic methods [10-16] as a function of actin monomer 
concentration. From these and other in vitro studies using var- 
ious labeling techniques the following picture has emerged of 
filament growth kinetics in the presence of ATP (see fig. [D- CD 
Monomers are added to a growing filament end as ATP-actin. 
(ii) Rapidly, the ATP is then hydrolyzed to ADP and phos- 
phate (Pi), both remaining bound to the monomer host (ADP- 
Pi-actin) [10, 14, 17-22]. A rate of 0.3s -1 was reported in ref. 
22 in the presence of Mg, assuming random hydrolysis uninflu- 
enced by neighboring monomers, (iii) After a long delay Pi re- 
lease into solution occurs, generating ADP-actin [23-25]. Re- 
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ported release rates are in the range 0.002 — 0.006s -1 [23-26]. 

A typical filament in a growth rate experiment is thousands 
of monomer units in length and thus consists mainly of ADP- 
actin. Hence the picture which emerges is of a long ADP-actin 
filament with a complex 3-state cap region at the filament end 
[5] (see fig. 0. A major goal of this report is to establish the 
composition and kinetics of the cap, and how these determine 
growth rates and measurable length fluctuations. This is im- 
portant in the context of cellular processes: the monomer com- 
position in actin filaments is thought to regulate actin-binding 
proteins in a timely and spatially organized way [2]. For exam- 
ple, it has been suggested that rates of branching generated by 
the Arp2/3 protein complex and/or debranching processes may 
depend on which of the 3 monomer species is involved, ATP- 
actin, ADP-Pi-actin or ADP-actin [7,26,27]. Pi release has 
been proposed to act as a timer for the action of the depolymer- 
izing/severing protein ADF/cofilin which preferentially attacks 
ADP-actin [2]. 

Our aim in this report is to establish theoretically the quan- 
titative implications of the currently held picture of actin poly- 
merization. Previous theoretical works addressed growth rates 
before the important process of Pi release was established [28- 
30]. To date there has been no theoretical analysis of single fil- 
ament non steady state growth rates rigorously accounting for 
the processes (i)-(iii) above. A recent theoretical work [31] has 
addressed steady state filament compositions. 

The cap has important consequences for the growth rate j 
as a function of ATP-actin concentration, c. Measured j(c) 
curves, such as those in fig. [51 are strikingly non-linear in 
the region near the concentration where growth rate vanishes 
[16, 32]. These become almost linear in excess Pi studies, 
where presumably the ADP-actin species is no longer involved 
[16]. The complexity of the cap structure and dynamics also 
underlies the values of the critical concentration c CI ;t at the 
fast-growing "barbed" end and slow-growing "pointed" end of 
the polar actin filament (c cr it denotes the concentration where 
mean growth rate at one end vanishes.) It is well known that 
in general these are different since detailed balance cannot be 
invoked for these non-equilibrium polymers [30]. Our work 
explores how these differences are related to cap structure. 

The major experimental focus has been mean growth rates, 
j(c). However, equally revealing are fluctuations about the 
mean whose measurement can expose features of the dynami- 
cal processes occurring at filament ends unavailable from j(c). 
These fluctuations are characterized by a "length diffusivity" 
D measuring the spread in filament lengths (see fig. \l\b)) 
similarly to simple one dimensional Fickian diffusion: after 
time t, the root mean square fluctuation in filament length is 
(2Dt) 1 / 2 about the mean value j(c)t. Using single filament 
microscopy, Fujiwara et al. [8] and Kuhn and Pollard [9] re- 
cently measured unexpectedly high values of this diffusivity 
near steady state conditions, D ~ 30mon 2 /s. This should 
be compared with what would be expected of an equilibrium 
polymerization involving the measured on/off rates of order 
lmon/s, which would lead to D ~ lmon 2 /s [8,30,33,34]. 
A number of possible explanations were proposed, (i) Fluctua- 
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Fi gure 1 : (a) Schematic of the 3-species cap at the barbed end of a long 
actin filament. Near the critical concentration we find a fluctuation-induced 
cap of N ca p Ri 25 monomers, with a short ATP-actin component, 7V^p P of 
order one. (b) Mean growth rate and fluctuations: in time t the average number 
of monomers added to a filament end is jt, with a spread of (2Di) 1 ^ 2 about 
this value. 



tions arise from "dynamic instability" due to stochastic cap loss 
episodes. This would be a far milder version of the "catastro- 
phes" in microtubule polymerization [8,35]. (ii) Filament poly- 
merization proceeds by addition and subtraction of oligomeric 
actin segment [8, 35]; this would constitute a radical departure 
from the accepted picture of filament growth kinetics involv- 
ing single monomer addition events, (iii) Growth involves ex- 
tra stochastic events such as short pauses possibly originating 
in filament-surface attachments [9]. (iv) Enhanced fluctuations 
results from an artifact due to monomer labeling [36]. (v) Ex- 
perimental error in filament length measurements [9]. A major 
focus of this report is to calculate the concentration-dependent 
length diffusivity, D(c), assuming that the standard monomer- 
by-monomer addition picture is valid. We will see that this 
leads to large D values below c CI ; t ; just above the critical con- 
centration fluctuations are enhanced, though much less than the 
experimental values. 

We consider the initial condition where long pre-formed 
ADP-actin seeds are exposed initially to a buffer of fixed actin 
concentration c and excess ATP. Thus for a given c value, a 
filament consists of a very long ADP-actin core at the end 
of which lies a complex steady state (but fluctuating) ATP- 
actin/ADP-Pi-actin cap. Our analysis emphasizes the barbed 
end, the pointed end assumed blocked. Our results apply to 
very dilute filaments where only ATP-actin is assumed to add 
to filaments since (i) free monomers bind ATP more strongly 
than ADP [37], and (ii) depolymerized ADP-actin or ADP-Pi- 
actin has enough time to exchange its nucleotide for ATP be- 
fore repolymerization. An important issue is the nature of the 
ATP hydrolysis mechanism: the experiments of refs. 20, 21 
support a random mechanism, though others have suggested a 
cooperative vectorial mechanism occurring at the interface be- 
tween ADP-Pi-actin and ATP-actin with rate 13.6s -1 [19,28]. 
In this study, random hydrolysis is assumed throughout. 

Parameter Values and Mathematical Methods 

One of the major aims of this work is to identify qualitative, 
but experimentally measurable, features of the growth kinet- 



ics which are independent of the precise values of rate con- 
stants, since the latter depend on experimental conditions such 
as ionic strength [38] and the values themselves are often con- 
troversial. The parameter values we use are shown in table 1 
in which is the polymerization rate constant of ATP-actin, 

and v^, Up and Vp are the depolymerization rates of ATP- 
actin, ADP-actin, and ADP-Pi-actin, respectively. The rates of 
ATP-hydrolysis and Pi release (both assumed irreversible) are 
rn and rpj, respectively. In addition, we will explore the ef- 
fects of changing some of these parameter values. Since the 
monomer at the tip makes bonds with two nearest neighbors, 
each belonging to a different protofilament, one expects that 
rate constants may also depend on the state of neighbors. Here, 
however, we study the simplest "one body" model, assuming 
on/off rates depend only on the attaching/detaching species [6] 
and that hydrolysis and Pi release rates are uniform along the 
filament. The influence of "many body" effects will be briefly 
discussed below. 

To calculate filament growth kinetics and composition one 
is faced with the formidable task of obtaining the steady state 
probability distribution of all possible actin monomer sequences 
along the filament: there are 3 possible states per monomer, 
so for filaments of N units long this necessitates solving 3 N 
coupled equations. We have managed, however, to obtain a so- 
lution for the mean elongation rate j(c) by projecting the full 
system of 3 N equations onto a set of just 3 exact equations for 
the return probabilities ipj ,ifjf, and ipf. These are the proba- 
bilities that a given monomer which was polymerized at t = 
is again at the tip at time t as ATP-actin, ADP-Pi-actin, or ADP- 
actin, respectively. 

The outline of our method is as follows. For j < the 
growth rate is related to the return probabilities by j = v^,p COIC , 
where p corc = 1 — f °° dt(ipf + tfjf + tpf) is the probabil- 
ity of exposure of the ADP-actin core at the tip. For j > 0, 
the relation is j = k^c — dtF t where F t = ipjv^, + 

ipfvp + ip^v^, is the mean depolymerization rate at time t of a 
monomer which added to the tip at t = 0. The integral of F t is 
the total depolymerization rate of added monomers. In the sup- 
porting material we present the dynamical equations obeyed by 
the return probabilities, from which we obtained a closed re- 
cursion relation for the Laplace transform of F t , namely 
This relates Je to fE+r a and /_E+ rpi . With boundary condi- 
tion Je — > as E — » oo, we started from large E values 
and evolved this equation numerically towards E = to ob- 
tain fo = J °° dtFt. Given /o, the time integrals of the return 
probabilities were then directly obtained from the dynamical 
equations and j was thereby determined. 

The above analytically based method does not generate cap 
sizes and length diffusivities. In order to calculate these quan- 
tities and also to test the validity of the analytical method we 
have simulated the stochastic tip dynamics employing the ki- 
netic Monte Carlo (MC) method known as the BKL [39] or 
Gillespie [40] algorithm, to evolve the state of a filament tip 
in time and to calculate its mean growth rate. Each step of the 
algorithm entails updating time by an amount depending on the 
rate and number of possible future events, namely polymerization/de- 
polymerization, hydrolysis, and Pi release. Excellent agree- 
ment is found between MC results and the numerical solutions 
of our closed equations for the growth rate (see inset of fig. 0*. 

Our analytical method is exact and avoids preaveraging, an 
approximation where the joint probability of a given filament 
nucleotide sequence is approximated as a product of probabil- 
ities for individual actin subunits. This neglects correlations 
between units. To assess the accuracy of this scheme, we com- 
pared our results for cap size and growth rate to those obtained 
using preaveraging (see supporting material for details). Preav- 
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Table 1 : Values of barbed end rate constants used in this work, appropriate 
for solutions of 50mM KC1 and ImM MgCl2. Units in s _1 unless otherwise 
indicated. ( a ) From ref. [6]. Assigned; at present there is no direct mea- 
surement of Up, but j(c) measurements with excess Pi [16] show the sum of 
the ADP-Pi-actin off rates at both ends together is a few times smaller than 
v~. < c ) From ref. [22]. < d > Fromrefs. [23-26]. 



eraging has been used in other theoretical studies of actin poly- 
merization such as ref. [31] to study steady state and ref. [32] 
to study growth rates. 

Cap Structure and the Importance of Fluctua- 
tions 

Using the parameters of table 1, in fig. [2 we present MC re- 
sults for (i) the total cap size, A cap , namely the mean total 
number of ATP-actin and ADP-Pi-actin subunits at the barbed 
end, as a function of concentration, and (ii) the number of ATP- 
actin cap subunits, N^ p . The figure shows that both caps be- 
come large for large concentrations. This is easy to understand. 
Consider for example the ATP cap: when polymerization rates 
exceed both the hydrolysis rate rn and the depolymerization 
rates, the interface between ADP-Pi-actin and ATP-actin fol- 
lows the growing tip with a lag of j(c)/rn monomers. Thus 
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Figure 2: Total cap length A/cap (c) and ATP-actin cap length N^ p (c) at 
barbed end. Parameters from table 1. Squares: MC results. Dashed lines: eq. 
ffl . Solid lines: preaveraging approximation. Vertical dashed line indicates 
c crit = 0.119/iM. 



mate dependence of cap length on concentration is 

,/[fc^(Ccrit-c)] , (c < C crit ) (2) 



A^cap = 



^cap^cap 



2D, 



cap / 



■j(c)/rn, N^=j(c)/r Pi (c» Ccrit ) • (1) 

Here the number of ADP-Pi subunits, JV^~ PPl , is found using 
similar reasoning as for A^p 13 . The validity of eq. Q for large 
concentrations is verified against MC data in fig. |2] 

The striking feature of fig. [2] is that the total cap remains 
long even below the critical concentration of the barbed end, 
being 25 units at c cr i t and remaining larger than unity down to 
c ~ c cr it/2. One might naively have guessed that below c cr it 
there would be no cap at all, since the filament is shrinking 
into its ADP core. (Indeed, the absence of a cap would also 
be suggested by eq. (0 if one were to extend its validity down 
to c cr it where j = 0.) This reasoning is however invalid be- 
cause it neglects fluctuations due to randomness of monomer 
addition/subtraction. 

To understand why fluctuations lead to long caps, consider 
the length changes of the cap only, excluding changes in the 
ADP-actin core length. Just below the critical concentration 
the tip of a typical long cap has a net shrinkage rate [33, 34], 
v cap (c). This is a weighted average of rates, summed over all 
possible states of the short ATP-actin segment on top of the 
long ADP-Pi-actin segment. Since v cap is a smooth function of 
c, it can be Taylor-expanded near the critical concentration and 
expressed as v cap = k^ s (c — c cr it) where k^ s is an effective 
on rate constant, different from fei. Now superposed on this 
average shrinkage the cap tip also performs a random walk in 
cap length space, described by a diffusivity -D C a P (c) [8,33,34], 
also an average over the states of the short ATP cap. (D cap is 
in fact the short-time diffusivity of the entire filament, see dis- 
cussion below.) For small times, diffusivity dominates and of 
order (2D cap t) 1 ^ 2 units add to or subtract from the cap. For 
times less than the cap turnover time < cap , this is much big- 
ger than the number of units wiped out by coherent shrinkage, Mean Growth Rate, j 

t 1 + I 'Via < • 'i n lirahtna ^ 1 o tVia hi-na »TfV»a»i ttio cnnnl/nna met 

^cap 



which indeed becomes large as c cr it is approached from below. 

In summary, even though on average below c cr it no ATP- 
actin monomers are being added to the tip, fluctuations in addi- 
tion/subtraction rates allow a cap to grow to length {2D cap t cap ) 1 / 2 
because the cap length diffusivity is dominant for times less 
than i C ap- Now since Pi release is very slow, for simplicity in 
deriving eq. (|2jl we assumed the release rate was zero, rp; = 0. 
However, the result of eq. (|2ji is valid even for a non-zero rp; 
except for concentrations so close to c cr it that the cap turnover 
time exceeds the Pi release time. In this inner region, diffusion 
is only able to grow the cap for a time of order rp, before Pi 
release intervenes. The maximum possible cap length, attained 
very close to c cr it, is thus 



Ng* « [2£cap(Ccrit)Api] 1/2 . (3) 

Eq. Q is valid until A cap reaches this bound. 

These arguments explain the origin of the long caps below 
Ccrit- To make quantitative comparison of eqs. 10 and Q to 
the numerics of fig. |2j the values of D cap and v cap must be 
determined. Now since for our parameter set and Vp have 
similar values (see table 1 ), an estimate can be obtained by con- 
sidering the special case where — Vp (identical ATP-actin 
and ADP-Pi-actin). This is a convenient case because D cap and 
Heap can be calculated exactly; the cap has just one monomer 

species, so = fef, and D C a P (c C rit) = (fc^c crit + v T )/2 
[8, 33, 34]. Using the values of table 1 in these expressions and 
in eq. gives N™!£ sw 26, of the same order as the numerics 

offig.m 

Finally, note that the preaveraging method shown in fig. 
13 is an excellent approximation in regions where fluctuations 
are unimportant (very large or very small c), producing almost 
identical results to MC. However, below c cr it it considerably 
underestimates cap lengths. This results from the preaveraged 
treatment of fluctuations. 



t. The cap lifetime < cap is the time when the shrinkage just 
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How is the behavior of the average rate of growth j(c) cor- 
related to cap structure and dynamics? The lowest curve of 
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Figure 3: Dependence of growth rate on concentration; Influence of 

(indicated in s -1 next to each curve). Other parameters as in table 1. MC 
and exact numerical solution results are indistinguishable. The spread in c cr i t 
values for the 3 curves is 5%. Left inset: blow up of critical region showing 
the agreement between MC (squares, error bars are standard deviation of mean) 
and numerical method (solid line). Right inset: Influence of many body effects; 
the value shown in s — 1 next to curves is depolymerization rate of ATP-actin 
next to ADP-actin, W^i D . 



fig. |3] shows numerical results for barbed end growth, using 
identical parameters to those of fig. |2] A noticeable feature is 
that the slopes are very different above and below the critical 
concentration of the barbed end. This directly reflects the cap 
structure just discussed, as follows. For c 3> c cr i t the ATP- 
actin segment is long and hides the remaining ADP-Pi-actin 
portion of the cap, so j w k^c — v^, has simple linear form 
and slope k^, approximately behaving as if ATP-actin were 
the only species involved. On the other hand for c < c cr it, 
the slope of j(c) is large because the cap length is changing 
rapidly as concentration increases (see fig. |2j. Filament length 
change is now generated by capless episodes, when the ADP- 
actin core is exposed and the filament shrinks with velocity 
(the steady state cap has fixed mean length and does not on av- 
erage contribute). Thus j = -v D p coro where p cmc « 1/-/V cap 
is the probability the cap length vanishes, assuming a broad 
distribution of cap lengths with mean N cap . Using eq. (|2j, this 

gives j » w E) fc ( ^ ff (ccrit - c)/(2D cap ) in the region where eq. 
(0 is valid. Since is large, this is a much larger slope than 
for concentrations above c cr ; t . 

The region very close to c cr ;t, where eq. takes over, 
is an interesting one: (1) Here the total cap becomes long, of 
length approximately N£%£, implying that ADP-actin is rarely 
exposed at the tip. This in turn implies that the depolymer- 
ization rate of ADP-actin will have only a small influence on 
the value of c cr it- This is verified in fig. [3] where we display 
j(c) curves for v D values ranging from 2.2 to 7.2s -1 . These 
changes produce only a very small shift in c cr ;t, even though 
j(c) changes significantly for c < c cr it- (2) The mean ATP-cap 
length is small (of order unity), and since the tip composition 
and cap length is constantly fluctuating, both ATP-actin and 
ADP-Pi-actin are frequently exposed at the tip. Thus we ex- 
pect a dependence of c cr it on the value of Vp . This is verified 
in fig. |4]where we display how the growth rate and c cr it change 
with the value of Vp . The magnitude of the shift is influenced 
by the assumed rate of ATP hydrolysis: if one uses, for exam- 
ple, a hydrolysis rate 10 times smaller, the change in growth 
remains substantial but is considerably reduced (see inset). 

Note also that preaveraging estimates the growth rate very 
accurately (see fig. 0). Even in the fluctuation-dominated re- 
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Figure 4: Growth rate: influence of the value of v~ (shown in s — 1 ). Other 
values as in table 1. Solid lines: numerical solutions and MC simulations (in- 
distinguishable). Dashed line: preaveraging approximation for v~ = 1.1s . 
Inset: same but with rjj = 0.03s — 1 . 



gion just below c cr it, where cap size is substantially underes- 
timated, it remains accurate though slightly less so than else- 
where. 

An important question is the effect of many body interac- 
tions between actin subunits, so far neglected in this report. 
We have found that the shape of the mean growth rate near 
and below the critical concentration is sensitive to these. As an 
example, fig. |3](inset) shows the dependence of j(c) on the de- 
polymerization rate of ATP-actin when its nearest neighbor is 
ADP-actin (w^| D ), with all other rates as in table 1. Other types 

of many body interactions can lead also to shifts in c cr it (not 
shown). Including many body interactions rapidly increases 
the number of rate constants. Since these are unknown and 
presumably hard to measure, this limits the uniqueness with 
which growth rate curves can be modeled near c cr it- We stress, 
however, that the central qualitative conclusions, namely the 
existence of a long cap at c cr i t and the associated change of 
slope of the growth rate, are general. An example of fitting ex- 
perimental j(c) curves with a one body model is shown in fig. 



Fluctuations in Growth Rate 

Turning now to fluctuations in growth rates, we find these be- 
have dramatically around the critical concentration reflecting 
a mild version of the dynamic instability exhibited by micro- 
tubules [30,41]. In the inset of fig. [6]we used MC to evaluate 

the length diffusivity, D{t) = (( L 2 ) - ( L f)/(2t), where L 
is the number of subunits added/subtracted after time t, starting 
from filaments with steady state caps at t = 0. For c = 0.15/iM 
(above c cr it) we find D is essentially independent of time. Its 
magnitude is of order lmon 2 /s, as would be expected for a 
growth process of identical subunits which add/subtract with 
rates of order Is" 1 [8,30,33,34]. However, for c = O.l^M (be- 
low c cr it) D is increasing with time, reaching a large asymp- 
totic value after several hundred seconds. Fig. |6|shows the 
dependence of D x on concentration; it exhibits a sharp peak 
below Ccrit and then drops rapidly. 

To understand the physics underlying this behavior, con- 
sider the simple model where ATP-actin and ADP-Pi-actin are 
identical (v^ — Vp) and Pi release very slow (rpi — > 0). 
Now D describes the random walk performed by the filament 
tip; if the tip makes a random forwards or backwards step of 
L monomer units every time interval T, then one can write 
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Figure 5 : Growth rate j (c) versus concentration from data taken from fig. 
1 of ref. 14 for simultaneous growth at both ends (in KC1 and Mg). Solid 
line: numerical results, barbed end (parameters from table 1), multiplied by a 
prefactor to fit data which lack absolute scale. Differences between numerical 
and experimental results may originate from the pointed end contribution, or 
possibly due to the experimental ionic conditions. 



D = L 2 /T. Just above the critical concentration, where on 
and off rates are approximately equal, the tip randomly adds 
or subtracts one ATP-actin (L = 1) in a mean time T = 
giving D = v^. Just below the critical concentra- 
tion, however, we know there is a long steady state cap. Since 
most filaments are capped, at short times D is determined by 
length changes of the cap and its value is thus close to the cap 
diffusivity, D cap . As time increases, more and more uncap- 
ping episodes occur, each episode now contributing to filament 
length change. Such events are correlated on the timescale of 
the cap lifetime, t cap w N 2 ap /v^. (we used D cap = for the 
simple model [33,34].) This explains why D(t) changes with 
time up to the cap lifetime (see fig. [6] inset). Thus to determine 
Z>oo, one must take T = t cap . Using a well known result from 
the theory of ID random walks [42], the number of uncapping 

events during the time t cap is approximately (Dcapicap) 1 ^ 2 ~ 
A^cap. Since the number of core monomers lost during each un- 
capping episode before a polymerizing monomer arrives is of 
order v^/v™, thus L = N cap v^/v^. This leads to a very dif- 
ferent expression for the diffusivity, Doc ~ (^d) 2 /' 1, t : Nereis 
a discontinuity in diffusivity at c cr i t of magnitude 

AD oc = v T (X 2 -1) , \ = v D /v T . (4) 

At the barbed end the instability parameter A w 5.1 and fluc- 
tuations at the critical concentration are very large, with a pro- 
nounced discontinuous drop in as one passes to higher c. 
A rigorous derivation of eq. @ is shown in the supporting 
material where in addition we obtain the full sawtooth curve 
shown in fig. |6j evidently, the simple model captures many 
features of the actual £>oo(c) profile. The effect of Pi release 
and ATP-actin/ ADP-Pi-actin differences is to shift c cr it and to 
smooth the sharp peak and shift it to somewhat below c cr it . 

How do the results of fig. [6] compare to the large fluc- 
tuations observed by Fujiwara et al. [8] and Kuhn and Pol- 
lard [9], and also suggested by the findings of ref. 43? Fig. [6] 
shows a peak value of D oo w 34mon 2 s~ 1 , dropping to Dqo w 
5mon 2 s _1 at c cr it- The experimentally reported value was ~ 
30mon 2 s _1 ; however, these measurements were performed at 
[8] or close [9] to a treadmilling steady state, i. e. at a concen- 
tration slightly above c cr it for the barbed end and well below 
that for the pointed end. At this concentration, fig. [6] shows 
a diffusivity less than 5mon 2 s _1 . Thus both theory and ex- 
periment exhibit large fluctuations near c cr it, but at different 
concentrations. Future experimental measurements of the full 
Doo (c) profile are needed to establish the relationship, if any, 
between these. 
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Figure 6: Long time length diffusion coefficient, Doo(c). Squares: MC 
results, parameters from table 1. Vertical dashed line indicates c cr i t . Solid 
line: prediction of simple model: v~ = v~ = 1.4s , rpi = 0, other values 
from table 1. Inset: Time-dependence of diffusivity at two concentrations. 



Our work leads also to the following prediction: since phos- 
phate will bind to ADP-actin and eliminate the effect of a large 
instability parameter, thus fluctuations and D at the barbed end 
will be suppressed in the presence of excess Pi. 

Discussion 

Pointed End j(c): Why is c cr it so Different? In this study 
we emphasized the barbed end, but our methods are also ap- 
plicable to the pointed end, provided the same mechanisms of 
uniform random hydrolysis and slow Pi release remain valid. 
Making this assumption, let us now discuss why c cr it (for ATP- 
actin) at the pointed end is almost six times the value at the 
barbed end [6]. Now an important issue is how different the 
ATP-actin and ADP-Pi-actin species are, in terms of on and off 
rate constants. That they are similar is suggested by the obser- 
vation that excess phosphate reduces the critical concentration 
in a pure ADP-actin polymerization to a value rather close to 
the barbed end c crit in ATP [16,44^-6]. But if indeed the 2 
species are similar, and the same basic mechanisms apply at 
the pointed end, this is inconsistent with the very different c cr ; t 
values. This inconsistency is due to the cap structure we have 
established here: the cap includes a long ADP-Pi segment es- 
sentially hiding the ADP-actin core which is thus rarely seen 
at the filament tip (see fig. 1(a)). For the barbed end (fig. 2) 
Acap ~ 25 at c cr it, and we find a large value for the pointed 
end at its c cr it, though smaller than the barbed end (data not 
shown). Thus ADP-actin on/off rates are almost irrelevant to 
c cr it (see fig. 3) and hence differences between ATP-actin and 
ADP-actin cannot account for the large c cr it differences. Thus 
the origin must be different ATP-actin/ADP-Pi-actin composi- 
tions at the pointed and barbed ends; since the ATP-actin seg- 
ment is short both species are regularly exposed at filament 
ends and substantially different c cr it values will result provided 
the 2 species have different rate constants. Were these identi- 
cal, Ccrit at both ends would be very similar, because the on/off 
rates at the filament ends would then be very close to the val- 
ues for an all ATP-actin filament; for such a filament, detailed 
balance dictates that the ratio of on/off rates at each end are 
identical [30]. (However, in apparent contradiction to this con- 
clusion are the findings of ref. [6] where different on/off ratios 
were reported at each end, under conditions where long ATP- 
actin caps are expected. A conceivable explanation is possibil- 
ity (ii), see below.) Many body effects will further affect c cr it- 
We are driven to two possibilities: (i) ATP-actin and ADP- 
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Pi-actin are substantially different, or (ii) different mechanisms 
operate during pointed end growth. Certain workers [47,48] 
have proposed possibility (i), based on the irreversibility of hy- 
drolysis [47] which suggests a large energetic change, possibly 
a structural change of the filament. Possibility (i) may in fact 
be consistent with the experiments of refs. 16,44^-6 which 
did not probe individual on/off rate constants of ADP-Pi-actin 
and which may have involved significant ADP-actin polymer- 
ization [45]. We are unaware of any crystallographic [49] or 
electron microscopic [50] experiments examining ATP/ADP- 
Pi differences for filamentous actin. 

If we adhere to the assumption that the growth mechanisms 
as previously outlined apply to both ends, we are then led to the 
following prediction: the values of c cr it for ATP-actin at both 
ends will be only weakly affected by the presence of excess 
Pi (provided ionic conditions are strictly unchanged). This is 
because the binding of Pi to ADP-actin segments is almost ir- 
relevant since these are rarely exposed at the tip due to long 
caps at c cr it. Indeed, for the barbed end no significant shift has 
been observed in the presence of Pi [16,44^-6]. For the pointed 
end, however, a reduction of c cr ;t has been reported in the pres- 
ence of Pi and barbed end capping proteins [16,44^-6]. This 
cannot be explained within the present framework and suggests 
possibility (ii). Future experiment will hopefully settle this im- 
portant issue. 

Conclusions. In this work filament growth rates j(c) and 
their fluctuations, as measured by the diffusivity D(c), were 
calculated as functions of ATP-actin concentration c. To our 
knowledge, this is the first rigorous calculation of these quan- 
tities accounting for all known basic mechanisms. Pantaloni 
et al. [28,29] studied j(c) at the barbed end in a work before 
the mechanism of Pi release was discovered. Infinitely fast Pi 
release and vectorial hydrolysis were assumed. Given the data 
available at that time, in order to explain the sharp change in 
slope of j(c) at c cr ;t (see e.g. fig. 5), they further assumed 
(i) strong three body ATP-actin/ ADP-actin interactions which 
lead to stable short ATP-actin caps, and (ii) zero hydrolysis rate 
of the nucleotide bound to the terminal monomer. In our work, 
the origin of the sharp change in slope is precisely the fact that 
Pi release is slow, similarly to an earlier model of microtubule 
polymerization [51]. 

Recently, Bindschadler et al. [31] studied the composition 
of actin filaments accounting for all three actin species at steady 
state. We have examined the preaveraging approximation used 
in their work and showed that it leads to very accurate j(c) 
curves, but the cap lengths are underestimated below c cr ; t . 

Here we have addressed random ATP hydrolysis only. Fu- 
ture work is needed to analyze the implications of the vecto- 
rial hydrolysis suggested by refs. 19,28. We showed that for 
random hydrolysis j(c) is linear far above the critical concen- 
tration. Growth rate experiments for both ends together in the 
absence of KC1 have exhibited non-linearities up to c = 10/iM, 
far above the critical concentration of the barbed end which is 
1 /iM under these conditions [10, 1 1]. In refs. 10, 28 this obser- 
vation was attributed to vectorial hydrolysis at the barbed end 
while in ref. 6 this was assigned to the non-linear contribu- 
tion of the pointed end whose critical concentration is xs 
under the same conditions. 

Perhaps our most interesting finding is that the long time 
diffusivity has a large peak below the critical concentra- 
tion c cr it of the barbed end, followed by a sharp drop in a nar- 
row range above c cr j t . This conclusion is quite general and its 
origin is the smallness of the Pi release rate and the large value 
of the off rate of ADP-actin at the barbed end. Future mea- 
surements of length diffusivities over a range of concentrations 
promise to provide new information and insight on the funda- 
mentals of actin polymerization. 

This work was supported by the Petroleum Research Fund, grant 



33944-AC7, and NSF, grant CHE-00-91460. We thank Ikuko Fuji- 
wara, Jeffrey Kuhn and Thomas Pollard for stimulating discussions. 



References 

[1] Bray, D. (2000) Cell Movements. From Molecules to Motility, 2nd Ed. 

(Garland Publishing, New York). 
[2] Pollard, T. D & Borisy, G. G. (2003) Cell 112, 453^465. 
[3] Pelham, R. J & Chang, F. (2002) Nature 419, 82-86. 
[4] Young, M. E, Cooper, J. A, & Bridgman, P. C. (2004) /. Cell. Biol. 166, 

629-635. 

[5] Korn, E. D, Carlier, M.-F, & Pantaloni, D. (1987) Science 238, 638-643. 
[6] Pollard, T. D. (1986) J. Cell Biol. 103, 2747-2754. 
[7] Amann, K. J & Pollard, T. D. (2001) Proc. Natl. Acad. Sci. USA 98, 
15009-15013. 

[8] Fujiwara, I, Takahashi, S, Takaduma, H, Funatsu, T, & Ishiwata, S. 

(2002) Nature Cell Biology 4, 666-673. 
[9] Kuhn, J & Pollard, T. D. (2005) Biophys. J. 88, 1387-1402. 
[10] Carlier, M.-F, Pantaloni, D, & Korn, E. D. (1984) J. Biol. Chem. 259, 

9983-9986. 

[11] Carlier, M.-F, Pantaloni, D, & Korn, E. D. (1985) J. Biol. Chem. 260, 
6565-6571. 

[12] Coue, M & Korn, E. D. (1986) / Biol. Chem. 261, 1588-1593. 

[13] Carlier, M.-F, Criquet, P, Pantaloni, D, & Korn, E. D. (1986) J. Biol. 

Chem. 261, 2041-2050. 
[14] Carlier, M.-F, Pantaloni, D, & Korn, E. D. (1986) J. Biol. Chem. 261, 

10785-10792. 

[15] Weber, A, Northrop, J, Bishop, M. F, Ferrone, F. A, & Mooseker, M. S. 

(1987) Biochemistry 26, 2537-2544. 
[16] Carlier, M.-F & Pantaloni, D. (1988) J. Biol. Chem. 263, 817-825. 
[17] Pardee, J. D & Spudich, J. A. (1982) J. Cell Biol. 93, 648-654. 
[18] Pollard, T. D & Weeds, A. G. (1984) FEBS Lett. 170, 94-98. 
[19] Carlier, M.-F, Pantaloni, D, & Korn, E. D. (1987) J. Biol. Chem. 262, 

3052-3059. 

[20] Ohm, T & Wegner, A. (1994) Biochim. Biophys. Acta 1208, 8-14. 

[21] Pieper, U & Wegner, A. (1996) Biochemistry 35, 4396-4402. 

[22] Blanchoin, L & Pollard, T. D. (2002) Biochemistry 41, 597-602. 

[23] Carlier, M.-F & Panataloni, D. (1986) Biochemistry 35, 7789-7792. 

[24] Carlier, M.-F. (1987) Biochemical and Biophysical Research Communi- 
cations 143, 1069-1075. 

[25] Melki, R, Fievez, S, & Carlier, M.-F. (1996) Biochemistry 35, 12038- 
12045. 

[26] Blanchoin, L, Pollard, T. D, & Mullins, R. D. (2000) Curr. Biol. 10, 
1273-1282. 

[27] Ichetovkin, I, Grant, W, & Condeelis, J. (2002) Curr. Biol. 12, 79-84. 
[28] Pantaloni, D, Hill, T. L, Carlier, M. F, & Korn, E. D. (1985) Proc. Natl. 

Acad. Sci. USA 82, 7207-7211. 
[29] Hill, T. L. (1986) Biophys. J. 49, 981-986. 

[30] Hill, T. L. (1987) Linear Aggregation Theory in Cell Biology . (Springer 
Verlag, New York). 

[31] Bindschadler, M, Osborn, E. A, Dewey, C. F. J, & McGrath, J. L. (2004) 

Biophys. J. 86, 2720-2739. 
[32] Keiser, T, Schiller, A, & Wegner, A. (1986) Biochemistry 25, 4899-4906. 
[33] O'Shaughnessy, B & Vavylonis, D. (2003) Phys. Rev. Lett. 90, 118301. 
[34] O'Shaughnessy, B & Vavylonis, D. (2003) Eur. Phys. J. E 12, 481^196. 
[35] Littlefield, R & Fowler, V. M. (2002) Nature Cell Biology 4, E209-E210. 
[36] Kudryasov, D. S, Phillips, M, & Reisler, E. (2004) Biophys. J. 87, 1 136- 

1145. 

[37] Kinosian, H. J, Selden, L. A, Estes, J. E, & Gershman, L. C. (1993) J. 

Biol. Chem. 268, 8683-8691. 
[38] Drenckhahn, D & Pollard, T. D. (1986) J. Biol. Chem. 261, 12754- 

12758. 

[39] Bortz, A. B, Kalos, M. H, & Lebowitz, J. L. (1975) Journal of Compu- 
tational Physics 17, 10-18. 

[40] Gillespie, D. T. (1977) / Phys. Chem. 81, 2340-2361. 

[41] Hill, T. L & Chen, Y.-D. (1984) Proc. Natl. Acad. Sci. USA 81, 5772- 
5776. 

[42] Fisher, M. E. (1984) J. Stat. Phys. 34, 667-729. 
[43] Brenner, S. L & Korn, E. D. (1983) J. Biol. Chem. 258, 5013-5020. 
[44] Rickard, J. E & Sheterline, P. (1986) J. Mol. Biol. 191, 273-280. 
[45] Wanger, M & Wegner, A. (1987) Biochimica et Biophysica Acta 915, 
105-113. 

[46] Weber, A, Pennise, C. R, & Fowler, V. M. (1999) J. Biol. Chem. 274, 
34637-34645. 

[47] Carlier, M.-F, Pantaloni, D, Evans, J. A, Lambooy, P. K, Korn, E. D, & 

Webb, M. R. (1988) FEBS Letters 235, 211-214. 
[48] Dickinson, R. B, Caro, L, & Purich, D. L. (2004) Biophys. J. 87, 2838- 

2854. 

[49] Otterbein, L. R, Graceffa, P, & Dominguez, R. (2001) Science 293, 708- 
711. 

[50] Belmont, L. D, Orlova, A, Drubin, D. G, & Egelman, E. H. (1999) Proc. 

Natl. Acad. Sci. USA 96, 29-34. 
[51] Hill, T. L & Carlier, M.-F. (1983) Proc. Natl. Acad. Sci. USA 80, 7234- 

7238. 



6 



SUPPORTING MATERIAL 

A. Numerical Method for Growth Rates. 

In this part of the supporting material we describe the numerical method we used in the main text to calculate the growth rate 
curves of figs. [3]and|4] Consider an ATP-actin monomer which polymerizes at the filament tip at time t = 0. We define the 
return probabilities ipj, ipf, and ijjf to be the probability that this monomer is once again at the tip after time t as ATP-actin, 
ADP-Pi-actin, or ADP-actin, respectively. The total depolymerization rate of this monomer at time t is 

F t = v T ^>J + v P $ + v D ^f . (5) 
The dynamical equations obeyed by the return probabilities are 

d f* 

ipj = -(k+c + + r H )^ T + k+c / dtty?*t-t'e _rH(t-< > , 



dt 
d 
df 



o 



= -(K c + v p + rpi)f^ + mi'J + k+c I dt'^Ft^e-^-^ 







dt'lpTrFt-t r -^— ( e-raC*-*') - e -™(f-*'A 

^Vf - -(4c + v D )V t D + rpi< + k+c J dt/1$F t -t> + k+cj^ dt' F t - V (l - e -^t*-*')) 

+ k+c f dt'^F t ^(l Ig—e-n^-O + rH e -rMt-n) . (6) 

Jo V r P i - r H rpi - r H / 

Here the non-integral terms on the right hand sides represent change of tip status due to polymerization, depolymerization, 
hydrolysis, and phosphate release events at time t. The integral terms represent rates of reappearance of the monomer at the 
tip, weighted by factors accounting for the probability of hydrolysis or phosphate release during the time interval since the last 
appearance at the tip. For example, the first term on the right hand side of the first equation represents the rate of change of 
the probability of finding the ATP-actin monomer at the tip due to (i) polymerization of another monomer on top of it, (ii) 
depolymerization of the monomer itself, or (iii) hydrolysis of the ATP nucleotide bound to the monomer at the tip. The integral 
term on the right hand side represents reappearance events of the ATP-actin unit at the tip given that it got buried inside the 
filament due to a polymerization event at time t', an event which occurred with rate k+>c. Factor F represents the rate of 
reappearance of the buried monomer at the tip due to depolymerization of all the monomers which were added on top of it. The 

factor e _rH( - t_ * ) is the probability that the ATP-actin monomer in question is not hydrolyzed while being buried. 
Now the filament growth rate is given by 

v D [l-f Q °°dt(^+tf+^)] (j<0) 

(7) 

k+c - / °° dtF t (j > 0) 

Carrying out a Laplace transformation, t — ► E, Ft — * /e, and ip t — > ^ e one has from eq. (0 

r t£[i a<o) 

3 = { (8) 

{ k+c - f (j > 0) 

while from eq. one obtains 

*1 = 1/ (E + v T +r H + k+c{l- f E+rH )) , 

= r H + k+c r H (f E +r H ~ /g+r Pi )/(r Pi - r H ) ipT 

E E + Vp + r Pi + k^c(l - fs+rpd El 

^ D _ (rpi + k+,c(f E - fE+r Pi )) g| + k+,c(f E - rpif E +r H /{rpi - r H ) + mfE+r Pi / {rpi - r H )) 

E ~ E + v D +k+c(l- f E ) 

Eliminating all ^ in the Laplace transformation of eq. (|5} after using eq. l|9) one obtains the following recursive relationship 
involving the function / alone: 

Se = T^[fB+mi /s+rpj , (10) 

where 



v\f f ]- ~ bl + V 6 " ~ 46260 nn 

l<-[jE+r H , JE+rpil — ■ I 11 ) 

Z09 
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Here the symbols bo, b\, and b 2 are functions of E, fE+r H and f E + rpi as follows: 

60 = A aa E 2 + Ao.jE + Aq.q , 

bi = A h3 E 3 + A h2 E 2 + A ltl E + A lfi , 

b 2 = A 2a E 2 + A 2 . 1 E + A 2 . a , (12) 

where 

A ,2 = -{ra -r Pi )v^k^c, 

Ao,i = (Vd^h - r Pi v^ + r H (-w P + v T )){k^c) 2 f E+rri + (Wp - v u r Pi )(k^c) 2 f E+ril 

- (th - r Pi )fc^c(r H Wp + % + w p + rpi + 2fc£ c)) 7 
A),0 = -v D (r H - r Pi ){k+cf f E+rF J E+ru 

- ( r a{v P - v T ) + r Pi v~)k+c + v D ((r H ) 2 - r Pi v T + r H (-r Pi + u~ + k+c))(k+c) 2 f E+rri 
+ (r H Upfc+c + WD(r H (up + r Pi ) - r Pi (v p + r Pi + fc+c)))(fc+c) 2 / B+rH 

- (r H - rp i )fc^c(fc^c(rHi'p + v^(vp + r Pi + k^cj) + v^(rn(vp + r Pi ) + v^(v P + r Pi + fcfc))) , 
4l,3 = r H - r Pi , 

4l,2 = -(m - r Pi )k+c(f E+rpi + f E+rii ) + (m - r Pi )(v D + r H + Vp + r Pi + t>~ + 3fefc) , 
= (m - ?'Pi)(fc^c) 2 /£; +r p i / £ ; +rH - (r H - r Pi )(u D + r H + w T + 2fc^c)fc r c/, E+rpi 

- (rH - r Pi )(t> D + Up + r Pi + 2fc^c)fc^c/ £+rH 

+ (?*H — r P i)(vpV^ + rpiUrp + 2vpk^c + 2r P \k^,c + 3v^k^c + 3(fcf c) 2 

+ v d ( r H + V P + r Pi + v t + c ) + r H(fp + rpi + 2fc^c)) , 
■^1,0 = (rH ~ r Pi ){v D + k+c)(k+c) 2 f E+rpi f E+rH 
+ (- u D( r H - r Pi v T + r H (-r Pi + v T + k^c)) 

+k^c(~ r H + r H (fp + r Pi - 2v^ - k^c) + r P i(2v^ + k^c)))k^cf E+rpi 
+ ( v d (- r n( v P + r Pi) + r P i(v p + r Pi + fc+c)) 

+k^c(r Pi (vp + r Pi + fc^c) - r H (2«p + r Pi + k^c)))k^cf E+ra 
+ (rH - rpi)(v D (r H (vp + rpi) + %(up + r Pi + fc+c)) 

+fc^c(rH(2wp + r Pi + fc^c) + (vp + r Pi + k^c)(2v^ + fcjc))) , 
^2,2 = r Pi - r H , 

^2,1 = (rH - r Pi )k^c(f E+rpi + f E+rii ) - (r H - r Pi )(r H + v p + r Pi + v T + 2k^c) , 
A 2 ,o = -(r H - rp i )(fc^c) 2 / B+rpi /_ E+ . rH + (r H - r Pi )(r H + + k^c)k^cf E+rpi 

+ (r H - r Pi ) (tip + r Pi + fc^c)fc^c/ B+rH - (r H - r Pi ) (tip + r Pi + fc^c)(r H + u T + fc^c) . (13) 

We remark that eq. ill 11 is the solution of a quadratic equation; which of the two solutions of the quadratic is meaningful is 
easily checked by demanding / < 1 in the E — > 00 limit. 

Now for any given monomer concentration c, with the boundary condition f E — > as E — > 00, we started from a large 
enough E value and evolved eq. (II 01 towards E = to obtain fo, f rpi , and /. rH . Substituting these values in eq. l|9} we 
further obtained "Pq , and 'J'g 3 . Thus, given /q, , and "Pq we evaluated j(c) using eq. (|8jl. It was shown that this 
method converges to a unique solution provided one starts the evolution from large enough E, retaining a sufficient number of 
significant digits. 

B. Preaveraging Approximation: Calculation of Growth Rates. 

In this part of the supporting material we present the method we used to calculate cap sizes and growth rates based on a 
preaveraging approximation. As discussed in the main text, compared to the exact calculations, this method gives different 
results for the cap size below c cr ;t, but provides very good approximations to growth rates. We denote <j^{n), 4> p (n) and 
</>d (n) the probability that the n th monomer away from the tip binds ATP, ADP-Pi or ADP, respectively. Consider first the tip, 
i.e. n = 1. One has [31] 

j t MV = *t c [M1) + Ml)] + b P <Ml) + v b MWt(2) - t-T [<M2) + 0d(2)]0 t (1) - raMl) ■ ( 14 ) 

The first term on the right hand side represents change of tip status into ATP-actin due to polymerization of ATP-actin at an ADP- 
Pi-actin or ADP-actin tip. The second term represents change of tip status into ATP-actin due to (i) depolymerization of ADP- 
Pi-actin or ADP-actin at n = 1, and (ii) exposure of ATP-actin, previously buried at position n = 2. Within the preaveraging 
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approximation, the joint probability of ADP at n = 1 and simultaneously ATP at n = 2, for example, is approximated as a 
product of probabilities: 0dt(1, 2) « </>d(1) </>t(2) in eq. (I14t . Similarly, the third term on the right hand side represents 
depolymerization of ATP-actin while the last term is change due to hydrolysis. 
For ADP-Pi-actin one has similarly 

= [v T <Ml) + VdMWpP) - {% [<M 2 ) + <M 2 )] + fcfc + ^Pi} Ml) + rH^r(l) • (15) 
The analogous equations for n > 1 are [31] 

1) - 0r(n)] + [%^t(1) + «d0d(1) + «p<Ml)][<M™ + 1) ~ 0t(»)] - r H <M™) 



1) - P (n)] + [u t 0t(1) + "dM 1 ) + % 0p(l)][0p(n + 1) - 0p(ti)] 

+ r H 0T("-) - r Pi P (n) (n > 1) (16) 

The rate of change of the ADP-actin probabilities are determined from (fc{nY + 6p(n) + 4>T>(n) = 1. Starting from an arbitrary 
nucleotide profile and using long filaments, we evolved numerically eqs. il4i-ll6i for a long enough time to allow the profile 
to reach its steady state (note that an analytical solution is also possible since j!6l> is linear, except for the tip terms [31]). The 
growth rate and the cap size were calculated from 

oo oo 

N cap = J2Mn) + Mn) , N^ F =J2Mn), 3 = k+c - v^M^) ~ v p ~ «d&>(1) • (17) 

n— 1 n— 1 



a 
It 



— (j)p(n) = fc^c[0 P (n 



C. Analytical Calculation of Length Diffusivity as a Function of Concentration. 

Here we prove that £>oo(c) in fig. [6]has a sawtooth shape in the special case Vp = and rpi — * 0. Consider first shrinking 
barbed ends, c < c cr j t . For long times, t 3> < cap , apart from fluctuations in the size of the steady state ATP/ADP-Pi cap, 
fluctuations in tip displacement are due to fluctuations in how far the ADP core has shrunk. Let us call u(r\t) the probability 
that in time t the filament is uncapped for a total time r. The probability that N monomers have been lost during this time is 

/>OC 

p(N,t)= dTii(T\t)T{N,v D T) , V(N,x)=x N e- x /N\, (18) 

Jo 

where the Poisson distribution, V, describes the probability distribution of the number of depolymerized ADP-actin monomers 
during the total uncapped period. Let us evaluate u by noting that the average value of r and its second moment are given by 



(T)=tS OQ , < r 2 ) = 2^ J df J* S(t" - t')dt" , 



(19) 



for long enough times. Here S(t) is a return probability, namely the probability that the tip is uncapped at time t, given that it was 
uncapped at t = 0, and S(t) — > Soo for t — > oo. To prove eq. ( |19> . define £(f'|f ) to be a random variable which is zero (unity) 

when the tip is capped (uncapped) at time t", given it was uncapped at f . One has ( r 2 ) = J ' f Q dt'dt" ( £(£'|0)£(t"|0) ) = 

2 /„' df f* dt" ( f (f |0) } ( £(f |f ) )■ Noting that S(t) = { £(t|0) ) one recovers eq. JT5). 
Now the exact result for S is [30] 



oc 



S{t) = e- {k ^ c+v ^ )t [h{2tx)+y~ 1/2 h{2tx) + {l-y)^y- j/2 I J {2tx)], y = k+c/v^ , x = {k+cv^) 1/2 , (20) 

where Ij are modified Bessel functions. Using Soo = 1 — c/c cr it one obtains for long times ( r ) = t(l — c/c cr ;t) and 
(t 2 ) c — (t/v^) where ( ) c denotes second cumulant. Thus relative fluctuations in r become small for long times and u 
becomes Gaussian, u(r\t) const. exp[— (t — ( r )) 2 /(2 ( r 2 ) )]. Substituting u in p(N, t) of eq. Jl 8I > and performing the 
integration one obtains 

D (c) = { ^d/ 2 )! 1 + Od/ w t - Wccrit] (c<c crit ) 
°° v j \ {k+c + v^)/2 (c > c crit ) 

which is the sawtooth curve plotted in fig. [6] Notice that decreases for smaller concentrations since at c = one must 
recover the Poissonian fluctuations of a pure depolymerization process for which Dqq = v^/2. In eq. (12 1 1 . the c > c cr it 
expression represents the fluctuations of a polymerization process of identical subunits [34] (since in the limit considered here 
the cap is never lost above c cr i t ). 
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